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Abstract 

A study has been performed that investigates parameter 
measurement effects on calculated in-flight thrust for the 
General Electric F404-GE-400 afterburning turbofan engine 
which powered the X-29A forward-swept wing research air- 
craft. Net-thrust uncertainty and influence coefficients were 
calculated and are presented. Six flight conditions were an- 
alyzed at five engine power settings each. Results were ob- 
tained using the mass flow-temperature and area-pressure 
thrust calculation methods, both based on the commonly 
used gas generator technique. Thrust uncertainty was deter- 
mined using a common procedure based on the use of mea- 
surement uncertainty and influence coefficients. The effects 
of data nonlinearity on the uncertainty calculation procedure 
were studied and results are presented. The advantages and 
disadvantages of using this particular uncertainty procedure 
are discussed. A brief description of the thrust-calculation 
technique along with the uncertainty calculation procedure 
is included. 


Nomenclature 


A 

cross-sectional area, in 2 

AE 

effective cross-sectional area, in 2 

ALT 

altitude, ft 

AP 

area pressure 

CFG 

gross thrust coefficient 

Ci 

influence coefficient 

E 

error in thrust due to parameter-measurement 
uncertainty 

ECU 

electronic control unit 

FG 

gross thrust, lbf 


* Aerospace engineer. Member AIAA. 

Copyright © 1989 by the American Institute of Aeronautics and Astro- 
nautics, Inc. No copyright is asserted in the United Slates under Title 17, 
U.S. Code. The U.S. Government has a royalty-free license to exercise all 
rights under the copyright claimed herein for Governmental purposes. All 
other rights are reserved by the copyright owner. 


FGI 

FHV 

FN 

FNAP 

FNWT 

FR 

FVG 

9 

H 

HPC 

IFT 

LPT 

LVDT 

M 

N 1 

N2 

PLA 

PT 

P 

R 

RSS 

T 

TT 

U 

UF 

UFNAP 

UFNWT 


ideal gross thrust, lbf 
fuel heating value, Btu/lbm 
net thrust, lbf 

net thrust, area-pressure calculation method, 
lbf 

net thrust, mass flow-temperature calculation 
method, lbf 

ram drag, lbf 

fan variable guide vane angle, deg 

gravity constant, 32.17 ft-lbm/lbf-s 2 

total enthalpy, Btu 

high pressure compressor 

in-flight thrust 

low pressure turbine 

linear-variable differential transformer 

Mach number 

fan rotor speed, rpm 

high pressure compressor rotor speed, rpm 

power lever angle, deg 

total pressure, lbf/in. 2 

static pressure, lbf/in. 2 

gas constant, lbf-ft/lbm-°R 

root-sum-square 

static temperature, °R 

total temperature, °R 

measurement uncertainty 

root-sum-square thrust uncertainty 

uncertainty in net thrust, area- pres sure 
method 

uncertainty in net thrust, mass flow- 
temperature method 


V velocity, ft/s 

VI ideal velocity, ft/s 

W mass flow, lbm/s 

WFAB afterburner fuel flow, lbm/s 

WFE engine core fuel flow, lbm/s 

WT mass flow temperature 

7 specific heat ratio, air 

F404 engine station identification numbers: 

0 free stream 

1 engine inlet 

2 fan inlet 

25 high pressure compressor inlet 

3 high pressure compressor discharge 

4 combustor discharge 

5 low pressure turbine discharge 

558 low pressure turbine discharge measuring 

plane 

6 afterburner inlet 

7 exhaust nozzle inlet 

8 exhaust nozzle throat 

9 exhaust nozzle discharge 

Introduction 

Knowledge of the uncertainty of calculated in-flight 
thrust is important in understanding the accuracy of aircraft 
performance values, including vehicle drag. 1 Errors in en- 
gine and aircraft parameter measurements required for the 
calculation of thrust directly affect thrust uncertainty. The 
degree of influence that these errors have depends on the 
characteristics of the specific engine model and thrust calcu- 
lation method utilized. Other sources of uncertainty, includ- 
ing engine model error, affect thrust uncertainty, but these 
effects were not considered in this investigation. 

Several studies have been undertaken in the past two 
decades at NASA Ames-Dry den Flight Research Facility in- 
vestigating thrust calculation methods and their correspond- 
ing sensitivity to measurement errors. These studies were 
performed for various engines installed in several types of 
aircraft, including the XB-70 2 and the F-lll. 3 

The primary purpose of this investigation is to docu- 
ment the net-thrust uncertainty and influence of measure- 
ment error for the F404-GE-400 afterburning turbofan en- 
gine. Problems in the analysis procedure, including influ- 
ence data nonlinearity, were also investigated to better un- 
derstand the uncertainty methodology used. 

The engine analyzed was installed in the X-29A forward- 
swept wing research aircraft (ship 1) which was flight tested 
at NASA Ames-Dryden. The results of this study are ap- 
plicable to other aircraft powered by the F404 engine which 


utilize an instrumentation system similar to that used for this 
analysis. Previous work involving the F404 engine installed 
in the X-29A includes a limited study of thrust sensitivity to 
measurement changes. A limited thrust uncertainty investi- 
gation was also performed, but a less accurate instrumenta- 
tion system was used. 4 

An in-flight thrust (IFT) computer program was used to 
calculate thrust for the F404 engine using measured parame- 
ter values as input. This program was supplied by the engine 
manufacturer. 5 The program uses two variations of the clas- 
sical gas generator thrust-calculation technique: the mass 
flow-temperature (WT) method and the area-pressure (AP) 
method. The error effects of ten important measured input 
parameters on thrust were determined for both methods. In- 
let ram recovery, bleed air, and horsepower extraction were 
estimated and used in the thrust calculation. Their effects on 
thrust uncertainty were not included in this analysis since the 
errors in estimating these aircraft installation factors were 
unknown. For a similar reason, engine model error inherent 
to the IFT program was not considered. 

Actual measured test data was not available for input into 
the IFT program at the time the analysis was performed. 
Therefore, this data was estimated using an F404 engine 
specification computer simulation program. 6 Actual X-29A 
aircraft and engine instrumentation accuracy data, required 
in the thrust-uncertainty analysis, was used in the study. 

This report presents parameter measurement effects on 
calculated net thrust for the F404 engine which was installed 
in the X-29A. Six simulated flight conditions, ranging from 
10,000 ft to 40,000 ft in altitude and from 0.4 to 1.6 Mach 
number (A/), were studied at throttle power lever angle 
(PL A) settings from part-power (70°PL,4) to maximum af- 
terburner power (130°PL^4). Thrust uncertainty was deter- 
mined using a common procedure based on the use of mea- 
surement uncertainty and influence coefficients. These in- 
fluence coefficient values were also calculated and are pre- 
sented. The effects of data nonlinearity on the uncertainty 
calculation were also studied. Results are presented and the 
advantages and disadvantages of using this particular un- 
certainty procedure are discussed. In addition, a brief de- 
scription of the thrust calculation technique and uncertainty 
calculation procedure is included. 

Engine Description 

The F404-GE-400 engine is a 16,000-lbf thrust class, low 
bypass, twin spool turbofan with afterburner. The engine 
incorporates a three-stage fan and a seven-stage high pres- 
sure compressor, each driven by a single-stage turbine. The 
fan and high pressure compressor guide vanes utilize vari- 
able geometry. Bleed air extraction is provided at the sev- 
enth stage of the high pressure compressor. The combus- 
tor is a through-flow annular type utilizing atomizing fuel 
nozzles. The afterburner can be fully modulated from mini- 
mum to maximum augmentation and uses fan discharge air 
and an afterburner liner to maintain a low engine skin tem- 
perature. The hinged-flap, cam-linked exhaust nozzle is hy- 
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draulically actuated. An engine accessory gearbox is driven 
by the compressor rotor. This gearbox drives the lubrica- 
tion oil pump, the variable-exhaust nozzle power unit, the 
generator, and both the main and afterburner fuel pumps. 
A schematic view of the F404-GE-400 engine with station 
designations is shown in Fig. 1 . 

The engine-control system consists of the throttle, main 
fuel control, electronic control unit (ECU) and afterburner 
fuel control. Throttle (power lever) movement is mechani- 
cally transmitted to a power amplifier which positions the 
main fuel control. During flight, PLA ranges from 31° 
(flight idle) to 130° (full power with afterburner). Interme- 
diate power (full nonafterbuming) occurs at 87° PLA . 

At power settings below intermediate, engine inlet total 
temperature (TT1) and throttle movement control the high 
pressure compressor rotor speed (N 2 ) through the main fuel 
control. At intermediate power and above, fan rotor speed 
(N\) is controlled by the ECU as a function of TT1, while 
N 2 remains essentially constant. The ECU senses engine 
and aircraft parameters, computes schedules, and maintains 
engine limits. The afterburner fuel control schedules fuel 
flow to the pilot and main spraybars. 

A single F404 engine is mounted in the fuselage of the 
X-29A and utilizes two side-mounted, fixed geomctir inlets 
optimized for transonic performance. 

Instrumentation 

The engine location of measurements used in the thrust 
calculation are shown in Fig. 1. 

A resistance temperature device is used to determine TT\ 
while N 1 is measured using an eddy-current instrument uti- 
lizing magnetic pickup from the fan. Fan variable guide 
vane angle ( FVG ) is measured using a linear-variable dif- 
ferential transformer (LVDT). 

The low pressure turbine discharge total pressure mea- 
surement, PT55S, a critical parameter in thrust calculation, 
is obtained using four five-element total pressure rakes. The 
20 PT55S pressures are measured by a multi-port differ- 
ential transducer referenced to a highly accurate absolute 
transducer. The differential transducer is thermally con- 
trolled by a heater-insulation blanket to maintain a con- 
stant temperature at which the unit was calibrated. The fi- 
nal PT55S value is an average of the 20 PT55S measure- 
ments. Values outside a specified tolerance are omitted from 
the average. The uncertainty in the PT55Z measurement is 
a root-sum-square (RSS) of the differential PT5 58 uncer- 
tainty and the reference-pressure uncertainty. 

The nozzle throat area, A8, is also measured using an 
LVDT. Volumetric flow meters are used to measure engine 
core and afterburner fuel flows ( WFE and W FAB respec- 
tively). Fuel temperatures are measured in both the gas gen- 
erator and afterburner fuel lines to permit conversion of vol- 
umetric values to mass flow. Fuel heating value (FHV) is 
a laboratory determined quantity. 


The thrust calculation also requires the measurement 
of freestream altitude (ALT) and Mach number (M). 
These measurements are obtained through the aircraft air- 
data system. 

Table 1 presents the range and absolute uncertainty for 
each measured parameter used in the thrust calculation. The 
measurement-uncertainty values were either supplied by the 
engine and instrumentation manufacturers or determined 
through in-house calibration. This instrumentation was in- 
stalled in the X-29A during the aircraft’s performance flight- 
testing phase. 


Calculation Procedures 

In-Flight Thrust Calculation 

Thrust Calculation Methods. Thrust was calculated 
according to two variations of the commonly used gas gen- 
erator technique: the mass flow-temperature (WT) method 
and the area-pressure (AP) method. The gas generator 
methodology is based on classical momentum, energy, and 
continuity laws. The primary difference between the WT 
and AP methods is due to the manner in which nozzle mass 
flow is calculated. 

The simplest form of the equation for ideal gross thrust 
(FGI) is based on the time rate-of-change of mass momen- 
tum at the nozzle of the engine. This equation assumes com- 
plete exhaust expansion to ambient pressure and therefore 
neglects the thrust pressure term at the nozzle exit. Using 
F404 station identification, the equation takes the form 


W 8 

— VI 9 
9 


where W 8 is the nozzle throat mass flow rate and V 19 is 
the ideal nozzle exit velocity. 

V 19 is expanded further in the above equation to ob- 
tain the mass flow-temperature equation in ideal gross thrust 
form. A gross thrust coefficient (CFG) is used to ob- 
tain actual gross thrust; and ram drag (engine inlet airflow, 
W 1, multiplied by freestream velocity, VO) is subtracted 
to give net thrust. The result is the net thrust mass flow- 
temperature equation 
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The nozzle mass llow term (W 8) is further expanded 
in the above equation to give the net thrust area-pressure 
equation 


FNAP = 


|aF 8 PT8 7 




x 



(PT 8\ 

V pO ) 



w\ 

x CFG VO 


( 2 ) 


The CFG modifies the ideal gross thrust term of the equa- 
tions to account for incomplete expansion, two-dimensional 
expansion effects, nozzle friction, and other sources of 
thrust loss. The values for this coefficient are deter- 
mined empirically through engine altitude facility and 
ground-testing. 

Equations (1) and (2) illustrate the effects that nozzle 
and frecstrcam parameters have on the thrust calculation. 
The quantities in the equations are obtained through sev- 
eral measurements which are discussed in the following sec- 
tion. Complete derivation of both equations is presented in 
Refs. 3 and 7. 

The In-Flight Thrust Program. The thrust values 
used to obtain the results of this report were acquired us- 
ing an in-flight thrust (IFT) calculation program. 5 This com- 
puter program was developed by General Electric for the 
Navy’s F404-F- 18 Propulsion System Integration Program. 

The IFT program determines mass flow, pressure, and 
temperature at the exhaust nozzle inlet by modelling the en- 
gine as a gas generator. The gas generator procedure uses a 
combination of engine performance models, engine compo- 
nent ground-test data, and actual measured engine and air- 
craft flight data. The IFT program uses these models and 
data to generate the values necessary for thrust calculation. 
The use of actual flight data allows the program to adjust for 
engine-to-engine performance variations. 

A schematic representation of the actual IFT calculation 
procedure and data flow is shown in Fig. 2. The 10 parame- 
ter measurements listed in Table 1 are used as input into the 
IFT program and are labeled in italics in the figure. 

Flight-condition measurements ALT and M, along with 
the measured engine parameters ATI, TT1, and FVG, are 
used in the inlet model and airflow calculation. Utilizing 
these parameters, W\ is calculated. Altitude (ALT) is used 
to calculate freestream pressure (pO) based on standard day 
correlations. This value, along with W 1 and M, allows the 
program to estimate inlet ram recovery using empirically de- 
rived data. Engine inlet total pressure (PT 1) is then calcu- 


lated. Freestream pressure (pO) is also used directly in both 
equations (1) and (2). 

Next, using an energy balance between the turbine and 
compressor, along with N\,W\, and PT 1, the total temper- 
ature, enthalpy, and pressure of the flow leaving the com- 
pressor (XT3, H 3, and PT3 respectively) are calculated. 
Compressor discharge airflow (W3) is determined by sub- 
tracting estimated nominal bleed airflow from W\. Horse- 
power extraction, also estimated by the program, is used in 
the energy balance also. 

The measurements WFE, FIfV, and the compressor 
model output are then used to calculate the energy rise across 
the combustor and turbine. Total temperature and enthalpy 
at the afterburner inlet, TT6 and // 6, are then determined. 
Total pressure at this station, PT6, is equal to the PT55H 
measurement. Afterburner inlet mass flow (IV 6 ) is calcu- 
lated by adding W FE to W 1 and subtracting bleed airflow. 

If the afterburner is not operating, TT1, HI, and W1 re- 
main unchanged from station 6. Nozzle inlet total pressure 
( PT1) is obtained by subtracting afterburner frictional loss, 
based on ground-test data, from PT6. If the afterburner is 
in operation TT1, HI, and W1 are then determined using 
FHV and WFAB measured values and station 6 condi- 
tions. Pressure losses because of heating and friction are 
included to obtain PT1. 

Since the flow is isentropically compressed from the noz- 
zle inlet to the throat, PT8 and TT8 equal PT1 and TT1 
respectively. Nozzle throat total pressure (PT8) is required 
in both equations (1) and (2) while TT8 is used in equa- 
tion (1). Nozzle inlet mass flow (IV7 ) is adjusted for noz- 
zle leakage giving mass flow at the throat, W8, which is 
required in equation (1), the mass flow-temperature thrust 
equation. The measured nozzle throat area, A8, is adjusted 
using a flow coefficient to give effective throat area, AE8. 
This value is required as input for equation (2), the area- 
pressure equation. Freestream pressure (pO), from ALT, is 
used to determine nozzle leakage and flow coefficients. 

Ideal gross thrust (FGI) is then calculated according to 
both the WT and AP thrust-calculation methods. The re- 
sults are adjusted using the CFG, which is determined from 
empirical data and is based on nozzle operating conditions. 
Actual gross thrust (FG) is then calculated. Freestream ve- 
locity (VO), calculated using M and ALT, is multiplied 
with W 1 giving FR. By subtracting this term from FG, net 
thrust (FN) is obtained for both thrust-calculation methods. 

Certain aircraft installation effects, including inlet 
spillage and nozzle drag, were not included in the calcula- 
tion because they are independent of the net-thrust calcula- 
tion procedure. 

Thrust-Uncertainty Calculation 

The common procedure by which thrust uncertainty (U F) 
is calculated, is to perform a root-sum-square on the indi- 
vidual thrust errors (E) each due to the measurement unccr- 


4 


tainty of one parameter. 8 The thrust-uncertainty calculation 
equation therefore takes the form 

UF= (Ef + Bj + ...+ Ef 0 )' /2 (3) 

Equation (3) contains ten squared terms, one for each of 
the input parameters analyzed. 

The individual thrust-error values are determined ana- 
lytically by adding the associated measurement uncertainty 
(from Table 1) to the measurement-parameter value and in- 
putting the result into the IFT program along with the other 
required, but unmodified, parameter measurements. The 
calculated-thrust value is then compared to the value that 
results using all unmodified input, and E \ due to the single- 
measurement uncertainty, is determined. This procedure is 
repeated for each parameter. The E values must be recal- 
culated for a change in flight or engine condition or for a 
change in thrust-calculation method. 

Although precise, this thrust-uncertainty calculation pro- 
cedure can be time-consuming and cumbersome to use if 
multiple instrumentation systems, and therefore different 
measurement uncertainties, are analyzed. This is because of 
the need to rerun the IFT program for each change in a mea- 
surement uncertainty. Also, this procedure does not clearly 
differentiate between the influence that a parameter has on 
thrust and the resulting effect on thrust uncertainty. 

Because of these reasons, an alternate thrust-uncertainty 
calculation procedure is commonly used, based on the linear 
variation of sensitivity of thrust to changes, or uncertainty, in 
an input parameter measurement. The error in thrust ( E ) is 
estimated by multiplying the parameter influence coefficient 
(Ci) by the measurement uncertainty (U) of that parameter. 
The resulting equation takes the form 

UF= [(Cii x U]) 2 + (C» 2 x U 2 ) 2 + ... 

+ (Ciio x {7 10 ) 2 ] 1/2 (4) 

By definition, Ci is the slope of the data representing 
the change, or error, in thrust because of a change, or un- 
certainty, in a thrust-calculation parameter measurement. 
A large Ci value indicates a large parameter influence on 
thrust calculation. 

Taking advantage of the linear nature of the sensitivity 
data, each parameter Ci value was calculated by determin- 
ing the percent change in thrust due to a 1 -percent change 
in that parameter. This Ci calculation was accomplished 
by running the IFT program for each flight-engine condi- 
tion and varying the input of the parameter in question by 
this percentage. The resulting thrust value was compared 
to the baseline thrust value calculated using unmodified in- 
put, and a percent change in thrust was obtained. Because it 
is defined as a slope, the Ci, in percent form, was numeri- 
cally equal to the percent change in thrust since a 1 -percent 
measurement parameter change was used. By changing the 
flight condition, engine power setting, or thrust calculation 


method, a parameter’s Ci value also changed and had to 
be re-calculated. 

Since the Ci values were calculated in percent form, U 
values also had to be nondimensionalized for each flight 
and engine condition. This nondimensionalization was per- 
formed by dividing the absolute uncertainty value (from Ta- 
ble 1) by the actual measured value of that parameter. The 
calculated UF was then also in percent form. The percent 
UF value is relative to the baseline-thrust value. 

An advantage of using the thrust-uncertainty calcula- 
tion procedure, based on equation (4), is that once the Ci 
database is obtained, the IFT program is no longer needed. 
This is because changes in instrumentation only affect the 
parameter’s U value, not the Ci, which is inherent to the 
thrust-calculation method utilized. Also, availability of 
the Ci data allows thrust-uncertainty calculation to be per- 
formed despite unavailability of the IFT program. This 
thrust-uncertainty calculation procedure is commonly used 
and is the procedure on which the results in this report 
are based. 

Results and Discussions 

Thrust-uncertainty and influence-coefficient values were 
calculated for both the WT and AP method net-thrust results. 
Six simulated flight conditions were analyzed: 10,000 ft, 
M = 0.4 and 0.8; 30,000 ft, M = 0.9 and 1 .2; and 40,000 ft, 
M = 0.8 and 1.6. Each flight condition was analyzed 
at five engine PL A settings each: 70° (part-power), 87° 
(intermediate), 92.5° (minimum afterburner), 109° (mid- 
afterburner), and 130° (maximum afterburner). 

Flight data was not available for input into the IFT pro- 
gram at the time this analysis was performed. Because of 
the unavailability of flight data, an F404 engine specification 
program was used to estimate the required IFT program in- 
put measurement data. 6 Unlike the IFT program, the specifi- 
cation program is a complete simulator requiring only flight 
condition and engine power setting as input. It estimates op- 
erating parameter values throughout the entire engine. Like 
the IFT program, the specification program is based on ex- 
tensive engine ground and altitude cell testing. 

Calculated Thrust Uncertainty 

Figure 3 displays net-thrust uncertainty against PLA for 
all six flight conditions analyzed. These values were calcu- 
lated using the thrust-uncertainty procedure based on equa- 
tion (4). Figure 3(a) presents results using the WT thrust- 
calculation method Ci values and Fig. 3(b) displays results 
using AP method Ci values. 

Figure 3(a) shows that WT net-thrust uncertainty 
( UFNWT ) values generally decrease (i.e., improve) as 
PLA increases to intermediate power then increase to a 
peak at mid-afterburner and decrease again as maximum 
power is approached. Figure 3(b) shows that the AP net- 
thrust uncertainty ( UFNAP ) values generally improve 
steadily as PLA is increased. The lowest uncertainty values 
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normally occur at maximum power since many of the pa- 
rameter values are at their lowest measurement uncertainty 
(in percent) at full-power. 

All thrust-uncertainty values calculated are between 
1 percent and 11 percent of the baseline-thrust values. The 
average UFNWT value is 3.19 percent while the average 
UFNAP value, 6.11 percent, is nearly twice as large. The 
lowest UFNWT value calculated, 1.16 percent, occurs at 

10,000 ft, M = 0.4, and 130°PLA; the highest, 7.99 per- 
cent, occurs at 40,000 ft, M = 0.8, and 109 °PLA. The 
lowest UFNAP value, 3.00 percent, was also calculated at 

10.000 ft, M = 0.4, and 130 °PLA\ the highest, 10.41 per- 
cent, at 40,000 ft, U = 1 .6, and 70° PL A. 

At the 10,000 ft conditions, the uncertainty values change 
more gradually with PL A change for both thrust calcula- 
tion methods than at the higher altitudes and are generally 
smaller in magnitude. 

The data show that, at the same altitude, the higher Mach 
number UFNWT values are generally less than the val- 
ues at the lower M condition. This trend is reversed for the 
UFNAP data. 

Except for the 40,000 ft, M = 0.8, 109°PLA condition, 
UFNAP values are greater than the UFNWT values at 
the same engine operating condition. These results indi- 
cate that the mass flow- temperature method is more accu- 
rate than the area-pressure method in calculating thrust for 
the in-flight instrumentation sysLem analyzed. 

Calculated Influence Coefficients 

Figures 4(a) through 4(1) present influence coefficient 
data for each flight condition and thrust-calculation method. 
To compare the parameter Ci values relative to each other, 
the Ci data are presented against PLA in the form of area 
graphs, with each Ci value displayed relative to the shaded 
region below it, not relative to the PLA axis. Tables 2(a) 
through 2(0 present the Ci data in numerical form. Each 
table contains Ci data for a single-flight condition as well 
as the absolute net-thrust values calculated by the IFT pro- 
gram using both the WT and AP methods. 

The majority of the calculated Ci values are less than 

1.0 percent, except for changes in PT558 and AS in the 
AP method, which generally produce higher values in the 
1-pcrccnt to 3-percent range. In general, the Ci values cal- 
culated for the same parameter differ only slightly between 
thrust calculation methods at the same flight condition- 
PLA setting. However, PP558 and AS Ci values, us- 
ing the AP mediod, are generally much greater than the WT 
method PT55S and AS Ci results (on the order of 1 percent 
to 2 percent larger). Ci trends according to PLA setting 
vary considerably between different parameters. 

All data considered, FHV and PT55S have the largest 
average Ci values for the WT method. For the AP method, 
PT 558 and AS Ci values are largest on average. FHV di- 
rectly affects combustion temperature, and therefore TT8, 
and so it has a large influence on the WT method (eq. 1). 


Area-pressure (AP) PT 558 and AS Ci values are large be- 
cause of the predominance of PT 8 and AES in the FNAP 
algorithm (eq. 2). Fan rotor speed (N 1) is quite influential 
in the AP method but has minimal effect in the WT method 
due to sensitivity cancellation occurring because of the pres- 
ence of mass flow terms in both the FG and FR terms. Fan 
variable guide vane angle (FV G) has the smallest Ci values 
on average. 

The large influence of both PT55S and AS on the AP 
method indicates that these two parameters must be mea- 
sured w ith high precision to keep AP method-thrust uncer- 
tainty at a low value. 

Nonlinearity Analysis* The thrust-uncertainty calcula- 
tion method using Ci values will lose accuracy if nonlin- 
earity exists in the thrust change data, implying a varying 
Ci value. This is important to equation (4) since it assumes 
that the Ci does not vary with changing U values. However, 
most of the data exhibited at least slight nonlinearity. 

Some parameters, depending on the thrust-calculation 
method and flight and engine condition, display severe non- 
linearity. In certain instances, this has a noticeable effect 
on the thrust- uncertainty calculation when combined with a 
large percent measurement uncertainty according to equa- 
tion (4). 

Figures 5(a) and 5(b) present thrust change against pa- 
rameter change for two of the more severe cases of nonlin- 
earity observed. Figure 5(a) shows data for N 1 at 10,000 ft, 
M = 0.8, and 130°PLA. The curves show that both thrust 
calculation methods, WT in particular, become increasingly 
sensitive to negatively increasing changes in N 1 at this con- 
dition. Figure 5(b) displays data for FHV at 30,000 ft, 
Kf - 1.2, and 130°PLA. In this case, the AP method data 
display more linearity than the WT method data. In both 
of these cases, varying slope, and therefore nonconstant Ci, 
is obvious. 

In general, nonlinearity affects the WT method data to 
a greater extent than AP method data, with N 1 and FHV 
generally displaying the most nonlinearity. Mach number 
(A 4) exhibits more linearity than the other parameters, es- 
pecially at larger parameter changes. This is not surprising 
since the ram drag term is a linear function of M. Consistent 
data linearity is exhibited by PT 558 as is shown in Fig. 5(c) 
for the 30,000 ft, M = 0.9, and S1°PLA condition. Data for 
the -5 percent parameter change condition are not shown 
since the IFT program had difficulty running this point. 

To make sure that the constant Ci assumption gave ac- 
ceptable thrust-uncertainty results, thrust-uncertainty val- 
ues calculated using equation (4) were compared to val- 
ues calculated using equation (3) (the more precise, but 
time-consuming, procedure). The results are presented in 
Fig. 6. For each flight-engine condition, the differences be- 
tween the two thrust uncertainty values are plotted against 
values calculated using equation (4). As can be seen, the 
data show little difference between the procedures. For the 
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WT method, ihe absolute average discrepancy between the 
different thrust uncertainty calculation procedure values is 
0.11 percent; for the AP method, the average difference, 
0.04 percent, is even smaller. One point on the figure dis- 
plays a difference of approximately -0 .9 percent between 
the two procedures and occurs at 40,000 ft, M = 0.8, and 
92.5 °PLA for the WT method. The difference is mostly 
because of measurement effects of WFAB. The percent 
measurement uncertainty for this parameter (64.4 percent) 
is very large because of the low afterburner fuel flow rate at 
this condition. Because of this, even the slight nonlinearity 
that exists for WFAB causes the difference seen. 

The overall effects of nonlinearity are small. The re- 
sults validate the commonly used uncertainty calculation 
procedure based on equation (4) as a means for easily 
and accurately analyzing frequent changes to an instrumen- 
tation system. 

Measurement Uncertainty Contributions to 
Thrust Uncertainty 

As described in the Calculation Procedures section, the 
thrust uncertainty equation is the root-sum -square of sev- 
eral thrust errors, each due to the measurement uncertainty 
of a specific parameter. Each thrust error contributes to the 
thrust uncertainty value, with the contribution equaling the 
Ci value multiplied by the percent parameter uncertainty. 
Results show that a large portion of each root-sum-squared 
thrust uncertainty value can generally be attributed to just 
one or two of these individual parameter contributions. This 
is normally because of an exceptionally large Ci for that pa- 
rameter, or the parameter uncertainty is much larger than the 
others, or both. 

All of the thrust uncertainty contribution data was ana- 
lyzed, but because of the large quantity of data, only a rep- 
resentative example is presented. Parameter values, mea- 
surement uncertainty in percent, Ci values, and individ- 
ual parameter thrust uncertainly contributions in percent arc 
presented in Tables 3(a) and 3(b) for the 30,000 ft, M = 
0.9, 81°PLA and 30,000 ft, M = 1.2, 130 °PLA conditions 
respectively. Also shown are the root-sum-square thrust- 
uncertainty values. The WT method absolute net-thrust val- 
ues for these conditions, calculated by the 1FT program, are 
shown in the table subheadings. 

In the intermediate power case, Table 3(a), WFE has 
the largest affect on thrust uncertainty for the WT method 
due to its rather large measurement uncertainty, 5.23 per- 
cent, coupled with its fairly large Ci value, 0.46 percent. 
Nozzle throat area (>18) is the primary error source for the 
AP thrust-uncertainty value. Its measurement uncertainty, 
4.66 percent, is also large as is the Ci va'uc for this condi- 
tion, 1.61 percent. Even though the parameter uncertainty 
value for FVG is large, it has virtually no influence on thrust 
because the Ci value is very small (less than 1 x 10 4 per- 
cent) for this condition. In the maximum power case, Ta- 
ble 3(b), WFE is no longer the primary source for the WT 
thrust uncertainty value; WFAB and TT 1 are responsible 


for a large portion of the thrust uncertainty with the intro- 
duction of afterburner fuel flow and a higher Ci value for 
TT\ (0.44 percent compared to 0.27 percent in the inter- 
mediate case). For the AP thrust uncertainty, A8 is still the 
primary uncertainty source, but has less effect than in the in- 
termediate power case because of a lower nozzle area uncer- 
tainty, 2.57 percent, coupled with a slightly lower Ci value, 

1 .38 percent. 

Considering all of the data, A8 and PT 558 produce the 
largest contributions to AP method thrust-uncertainty val- 
ues on average. For WT method thrust-uncertainty values, 
the fuel flows, WFAB and WFE, produce the largest con- 
tributions. Engine inlet total temperature (TT1) contribu- 
tions are substantial for both methods also. Fan rotor speed 
(ATI) contributions are significant for the AP but not the 
WT thrust-uncertainty values. Despite the fact that their 
Ci values are significant, FHV and ALT contributions are 
small for both methods. This is because of their small mea- 
surement uncertainty (high measurement accuracy). The 
smallest contributions on average are produced by M, FV G, 
FHV, and ALT. 

Engine core fuel flow (WFE) and WFAB both have 
average WT method Ci values relative to other parameter 
results. However, the percent measurement uncertainty is 
quite large for both, especially at lower PLA settings (and 
hence lower fuel flow rates), and so their uncertainly contri- 
butions arc substantial for the WT method values. 

Concluding Remarks 

The measurement effects on the calculation of in-flight 
thrust for an F404-GE-400 afterburning turbofan engine 
were documented in this study. The mass flow-temperature 
(WT) and area-pressure (AP) thrust-calculation methods 
were used to calculate the thrust values used in the analy- 
sis. Six flight conditions throughout the flight envelope of 
the engine were analyzed at five engine power settings each. 

One of the primary purposes of this investigation was 
to document net-thrust uncertainty using an instrumenta- 
tion system flown on board the X-29A research aircraft. 
The analysis revealed that the average thrust uncertainty 
due to measurement uncertainty using the WT thrust cal- 
culation method was 3.19 percent while the average using 
the AP method, 6.11 percent, was nearly twice as large. 
All of the thrust-uncertainty values fell between 1 percent 
and 1 1 percent. For the WT method, the measurement ef- 
fects of engine core fuel flow, W FE and afterburner fuel 
flow, WFAB, made the largest contributions to thrust un- 
certainty. Nozzle throat area, .4 8, and turbine discharge to- 
tal pressure, PT558, contributed most to the AP method 
thrust uncertainly results. The lowest uncertainty values oc- 
curred at maximum power for both methods. Low altitude 
thrust-uncertainty values were generally less than higher al- 
titude values. 

To determine the influence of measurement error on the 
calculation of thrust, Ci values were calculated for all mca- 
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surcment parameters and for both thrust-calculation meth- 
ods. For the AP method, PT55S and AS Ci values av- 
eraged between 1 percent and 3 percent. All other param- 
eters generally had much lower values that averaged less 
than 1 percent for both methods. Fuel heating value (FHV) 
and PT55 8 had the largest Ci values, and therefore the 
most influence on net thrust, for the WT and AP meth- 
ods respectively. 

The effects of nonlinearity on influence coefficient and 
thrust-uncertainty calculation were investigated and found 
to be minimal except in a few isolated cases. 
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Table 1. Measurement ranges and uncertainly 
of thrust calculation parameters. 


Parameter 

Range 

Absolute measurement 
uncertainty, ± 

M 

0 to 2.0 

0.005 

ALT 

0 to 60,000 ft 

60.0 ft 

TTl 

400° to 860 °R 

8.6 °R 

N 1 

0 to 13,270 rpm 

132.7 rpm 

FVG 

0° to 55° 

i.r 

PT55H 

0 to 60 lbf/in. 2 

0.16 lbf/in. 2 

AH 

220 to 540 in. 2 

10.8 in. 2 

WFE 

0 to 12,000 lbm/hr 

240.0 lbm/hr 

WFAB 

0 to 30,000 lbm/hr 

600.0 lbm/hr 

FHV 

18,200 to 18,600 Btu/lbm 

93.0 Btu/lbm 
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Table 2. Influence coefficient and absolute thrust values. 
(a) 10,000 ft, U = 0.4 


Influence coefficient, percent 

WT method AP method 


Parameter 

70 

87 

PLA, deg 
92.5 

109 

130 

70 

87 

PLA, deg 
92.5 

109 

130 

M 

0.248 

0.193 

0.178 

0.151 

0.121 

0.248 

0.193 

0.184 

0.151 

0.120 

ALT 

0.201 

0.156 

0.152 

0.153 

0.158 

0.207 

0.169 

0.165 

0.160 

0.154 

TT\ 

0.680 

0.304 

0.315 

0.257 

0.440 

0.273 

0.129 

0.128 

0.123 

0.162 

N 1 

0.877 

0.189 

0.189 

0.144 

0.198 

0.525 

0.226 

0.223 

0.210 

0.223 

FVG 

0.048 

0.000 

0.000 

0.000 

0.000 

0.020 

0.000 

0.000 

0.000 

0.000 

PT 558 

0.437 

0.354 

0.339 

0.361 

0.436 

1.813 

1.704 

1.690 

1.662 

1.681 

A8 

0.159 

0.130 

0.125 

0.111 

0.088 

1.155 

1.082 

1.082 

1.066 

1.022 


0.433 

0.429 

0.373 

0.278 

0.099 

0.006 

0.005 

0.007 

0.015 

0.022 

W'PAB 

0.000 

0.000 

0.064 

0.271 

0.173 

0.000 

0.000 

0.000 

0.009 

0.014 

FHV 

0.426 

0.421 

0.434 

0.583 

0.329 

0.007 

0.003 

0.005 

0.026 

0.051 

Net thrust, lbf: 

4904 

7318 

7956 

9392 

11664 

4909 

7337 

7680 

9398 

11813 





(b) 10,000 ft, M 

= 0.8 





M 

0.550 

0.428 

0.388 

0.321 

0.244 

0.546 

0.427 

0.407 

0.324 


ALT 

0.238 

0.171 

0.164 

0.163 

0.160 

0.250 

0.188 

0.182 

0.168 

0.160 

TT 1 

0.602 

0.426 

0.435 

0.347 

0.549 

0.445 

0.380 

0.378 

0.344 

0.441 

7V1 

0.680 

0.451 

0.451 

0.336 

0.519 

0.937 

0.845 

0.830 

0.732 

0.816 

FVG 

0.078 

0.003 

0.003 

0.002 

0.004 

0.086 

0.004 

0.004 

0.003 

0.004 

PT558 

0.509 

0.425 

0.408 

0.410 

0.501 

2.238 

2.161 

2.141 

2.035 

2.123 

>18 

0.238 

0.097 

0.096 

0.106 

0.071 

1.452 

1.421 

1.392 

1.252 

1.215 

WFF 

0.537 

0.496 

0.422 

0.309 

0.103 

0.015 

0.043 

0.045 

0.046 

0.069 

WF>1B 

0.000 

0.000 

0.071 

0.2% 

0.210 

0.000 

0.000 

0.004 

0.023 

0.039 

FPV 

0.528 

0.487 

0.490 

0.636 

0.365 

0.016 

0.041 

0.047 

0.073 

0.137 

Net thrust, lbf: 

4724 

7928 

8750 

10564 

13892 

4760 

7957 

8334 

10467 

14299 





(c) 30,000 ft, M 

= 0.9 





M 

0.474 

0.387 

0.348 

0.281 


0.474 

0.387 

0.356 

0.277 


ALT 

0.567 

0.362 

0.362 

0.392 


0.595 

0.362 

0.365 

0.378 


TT 1 

0.735 

0.271 

0.278 

0.212 


0.806 

0.190 

0.184 

0.161 


7V1 

0.884 

0.086 

0.091 

0.061 


1.495 

0.347 

0.333 

0.287 


FVG 

0.020 

0.000 

0.000 

0.000 


0.028 

0.000 

0.000 

0.000 


PP558 

0.439 

0.411 

0.432 

0.408 


2.163 

1.960 

1.914 

1.829 

1.813 

>18 

0.044 

0.163 

0.125 

0.032 


1.552 

1.611 

1.526 

1.335 

1.234 

1VFF 

0.511 

0.461 

0.393 

0.295 


0.033 

0.068 

0.059 

0.046 

0.033 

1VF>1B 

0.000 

0.000 

0.072 

0.2% 

0.198 

0.000 

0.000 

0.009 

0.035 

0.025 

FHV 

0.502 

0.452 

0.462 

0.628 


0.031 

0.065 

0.066 

0.086 

0.076 

Net thrust, lbf: 

3395 

4547 

5061 

6261 


3399 

4553 

4944 

6359 

8145 
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Parameter 

70 

M 

0.526 

ALT 

0.329 

TTl 

0.320 

N 1 

0.134 

FVG 

0.000 

PT 558 

0.404 

A 8 

0.212 

WFE 

0.503 

WFAB 

0.000 

FIIV 

0.494 

Net thrust, lbf: 5475 


M 0.389 

ALT 0.727 

TT1 0.292 

Ml 0.107 

FVG 0.000 

PT5S8 0.441 

A8 0.024 

WFE 0.484 

WFAB 0.000 

FHV 0.476 

Net thrust, lbf: 2348 


M 0.679 

ALT 0.521 

TT 1 0.473 

Ml 0.281 

FVG 0.007 

PT558 0.553 

-48 0.297 

WFE 0.565 

WFAB 0.000 

FHV 0.554 

Net thrust, lbf: 4251 


Table 2. Concluded. 

(d) 30,000 ft, M = 1,2 

Influence coefficient, percent 
WT method 


P LA, deg 
92.5 109 


87 

0.526 

0.330 

0.320 

0.138 

0.000 

0.404 

0.210 

0.504 

0.000 

0.494 

5478 


0.344 

0.539 

0.258 

0.054 

0.000 

0.415 

0.127 

0.455 

0.000 

0.446 

2675 


0.679 

0.521 

0.473 

0.281 

0.007 

0.553 

0.297 

0.565 

0.000 

0.554 

4252 


0.466 

0.325 

0.328 

0.156 

0.000 

0.389 

0.178 

0.423 

0.074 

0.494 

6174 


0.312 

0.549 

0.250 

0.053 

0.000 

0.431 

0.084 

0.400 

0.079 

0.472 

2954 


0.593 

0.493 

0.477 

0.308 

0.007 

0.580 

0.233 

0.459 

0.085 

0.541 

4872 


0.371 

0.340 

0.250 

0.102 

0.000 

0.382 

0.090 

0.305 

0.298 

0.636 

7766 


0.459 

0.476 

0.313 

0.148 

0.004 

0.532 

0.163 

0.324 

0.326 

0.677 

6273 


130 

0.275 

0.359 

0.439 

0.237 

0.001 

0.432 

0.020 

0.108 

0.213 

0.371 

10477 


0.320 

0.431 

0.792 

0.759 

0.017 

0.481 

0.080 

0.102 

0.216 

0.368 

9012 


70 

0.544 

0.351 

0.402 

0.829 

0.001 

2.363 

1.933 

0.110 

0.000 

0.106 

5298 


0.391 

0.767 

0.265 

0.477 

0.000 

2.131 

1.457 

0.066 

0.000 

0.063 

2336 

0.687 

0.335 

1.774 

3.145 

0.052 

2.779 

2.407 

0.161 

0.000 

0.156 

4205 


(e) 40,000 ft, U = 0.8 


0.253 

0.198 

0.599 

0.618 

0.161 

0.352 

0.023 

0.051 

0.000 

0.000 

0.434 

0.429 

0.012 

0.047 

0.308 

0.124 

0.332 

0.232 

0.690 

0.399 

3639 

4654 


(f) 40,000 ft, M = 1.6 


AP method 
PL A, deg 
87 92.5 109 


0.541 

0.489 

0.372 

0.351 

0.341 

0.359 

0.401 

0.381 

0.336 

0.837 

0.785 

0.679 

0.001 

0.001 

0.001 

2.357 

2.256 

2.108 

1.926 

1.816 

1.550 

0.110 

0.098 

0.082 

0.000 

0.012 

0.051 

0.106 

0.109 

0.140 

5327 

5884 

7741 


0.347 

0.318 

0.250 

0.546 

0.558 

0.601 

0.179 

0.173 

0.154 

0.290 

0.280 

0.250 

0.000 

0.000 

0.000 

1.991 

1.945 

1.871 

1.519 

1.433 

1.253 

0.067 

0.058 

0.045 

0.000 

0.008 

0.031 

0.064 

0.064 

0.082 

2657 

2895 

3691 


0.685 

0.619 

0.453 

0.336 

0.360 

0.510 

1.771 

1.670 

1.375 

3.140 

2.936 

2.375 

0.052 

0.051 

0.038 

2.776 

2.685 

2.475 

2.403 

2.223 

1.891 

0.160 

0.158 

0.106 

0.000 

0.018 

0.082 

0.156 

0.174 

0.199 

4214 

4669 

6367 


130 

0.268 

0.399 

0.404 

0.644 

0.001 

2.118 

1.382 

0.093 

0.046 

0.169 

10740 


0.200 

0.606 

0.182 

0.240 

0.000 

1.836 

1.189 

0.038 

0.027 

0.082 

4598 


0.298 

0.518 

1.379 

2.256 

0.037 

2.432 

1.626 

0.081 

0.081 

0.213 

9683 
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Table 3. Parameter measurement values, measurement uncertainty, Ci values and 
thrust uncertainty contribution examples. 

(a) 30,000 ft, M = 0.9, 87 °PLA\ baseline FNWT = 4547 lbf 






Thrust uncertainty 


Baseline 

Parameter 

Influence Coefficient (Ci), 

contribution (Ci x U ), 

Parameter 

parameter 

uncertainty (17), 

percent 


± percent 


value 

± percent 

WT method AP method 

WT method 

AP method 

M 

0.9 

0.56 


0.387 

0.215 

0.215 

ALT 

30,000 ft 

0.20 

0.362 

0.362 

0.072 

0.072 

TT 1 

478.7 °R 

1.80 

0.271 

0.190 

0.486 

0.341 

N 1 

13,187.2 rpm 

1.01 

0.086 

0.347 

0.087 

0.349 

FVG 

0.00° 

* 

0.000 

0.000 

0.000 

0.000 

PT 558 

27.95 lbf/in. 2 

0.58 

0.411 

1.960 

0.238 

1.136 

AS 

231.9 in. 2 

4.66 

0.163 

1.611 

0.757 

7.502 

WFE 

4590.1 lbm/hr 

5.23 

0.461 

0.068 

2.408 

0.357 

WFAB 

0.00 lbm/hr 

0.00" 

0.000 

0.000 

0.000 

0.000 

FHV 

18,400 Btu/Ibm 

0.51 

0.452 

0.065 

0.229 

0.033 

Total RSS thrust uncertainty 

( UFNWT and UFNAP respectively): 

2.603 

7.615 

(b) 30,000 ft, M = 1,2, \30°PLA ; baseline FNWT - 10477 lbf 

M 

1.2 

0.42 

0.275 

0.268 

0.115 

0.112 

ALT 

30,000 ft 

0.20 

0.359 

0.399 

0.072 

0.080 

TT 1 

530.7 °R 

1.62 

0.439 

0.404 

0.710 

0.655 

Nl 

13,536.9 rpm 

0.98 

0.237 

0.644 

0.233 

0.631 

FVG 

0.24° 

458.33** *** 

0.001 

0.001 

0.273 

0.370 

PT 558 

34.27 lbf/in. 2 

0.47 

0.432 

2.118 

0.204 

1.001 

AS 

420.2 in. 2 

2.57 

0.020 

1.382 

0.053 

3.552 

WFE 

5652.2 lbm/hr 

4.25 

0.108 

0.093 

0.457 

0.394 

WFAB 

14107.1 lbm/hr 

4.25 

0.213 

0.046 

0.904 

0.198 

FHV 

18,400 Blu/lbm 

0.51 

0.371 

0.169 

0.188 

0.086 

Total RSS thrust uncertainty 

(i UFNWT and UFNAP respectively): 

1.325 

3.848 


♦The zero parameter value causes an infinite percent uncertainty value due to division by zero. The zero 
influence coefficient at this condition, however, causes zero thrust uncertainty contribution. 


**In this case, the percent uncertainty is zero since, with the valve shut off, it is known with certainty that 
there is no fuel flow. 

♦♦♦The high percent uncertainty for this parameter is due to the low measurement value. The very small 
influence coefficient at this condition for this parameter, however, causes a fairly small thrust uncertainty 
contribution. 
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Fig. 1 F404 engine station and measurement locations. 
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Fig. 2 F404 in-flight thrust calculation flow chart. 
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(b) AP method. 


Fig. 3 Net-thrust uncertainty as a function of PLA. 
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(a) WT method; 10,000 ft, M = 0.4. 



DU FHV ^ FVG 
I 'HA WFAB EH3 N1 

nun wfe E3 m 

m A8 ■■ ALT 
dl PT558 EED MACH 



70 80 90 100 110 120 130 


PLA, deg 

9259 

(b) AP method; 10,000 ft, M = 0.4. 

Fig. 4 Relative parameter Ci values. 
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(c) WT method; 10,000 ft, M = 0.8. 
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(d) AP method; 10,000 ft, M = 0.8. 
Fig. 4 Continued. 
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Fig. 4 Continued. 


120 


130 


16 






HE FHV E=3 FVG 
fZZl WFAB CD N1 

mm wfe ess m 

HH A8 MM ALT 
I I PT558 BSD MACH 





90 100 110 120 130 

PLA, deg 


9264 


(g) WT method; 30,000 ft, M = 1.2. 
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(h) AP method; 30,000 ft, M = 1.2. 


Fig. 4 Continued. 
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(j) AP method; 40,000 ft, M = 0.8. 
Fig. 4 Continued. 
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(k) WT method; 40,000 ft, U - 1.6. 
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(1) AP method; 40,000 ft, M = 1.6. 

Fig. 4 Concluded. 
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(a) Nl; 10,000 ft, M = 0.8, 130°PLA 
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(b) FHV\ 30,000 ft, M = 1.2, 130 °PLA. 

Fig. 5 Nonlinear data analysis examples. 
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(c) PT55 8; 30,000 ft, M = 0.9, &TPLA. 

Fig. 5 Concluded. 
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Fig. 6 Nonlinearity effects on equation (4) thrust uncertainty results. 
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